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Abstract. A detailed study of the absorbing molecular 
clouds towards the radio source PKS1413+135, at a red- 
shift z=0.247, is reported. Physical conditions (density, 
temperature, filling factors) are derived for the molecu- 
lar gas, in the frame work of two models: a homogeneous 
multicomponcnt model with equal filling factors, and a 
two-phase medium consisting of dense clumps embedded 
in a more diffuse component. It is shown that our absorp- 
tion data are consistent with the presence of a diffuse gas 
component, dominating the observed opacity, and a dense 
component, accounting for most of the mass. We also show 
that without knowledge of the small scale structure of the 
absorbing molecular gas, we can only derive lower limits 
to the column density. Given the very narrow absorption 
spectrum, the size of the overall absorbing cloud along the 
line of sight must be quite small, of the order of 1 pc. The 
variability of the absorption spectrum has been studied 
over a time range of more than 2 years. We find that the 
opacity ratio between two absorbing components in CO 
has varied by a factor 2.3 ± 0.3. The variations are inter- 
preted as a change of the line of sight due to structural 
changes in the background source. We discuss what can be 
derived from molecular absorption line observations con- 
cerning the invariance of physical constants, such as the 
mass of the molecules. We discuss molecular line data from 
redshifts 0.27-0.89, corresponding to look-back times of 
30% to 60% of the age of the universe. 

Key words: interstellar medium: molecules - interstel- 
lar medium: clouds - galaxies: BL Lac objects individual: 
PKS1413+135 - galaxies: ISM - galaxies: quasars: absorp- 
tion lines - radio continuum: galaxies 
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1. Introduction 

Studies of absorption lines towards high redshift QSOs 
give detailed information about the interstellar medium 
(ISM) in distant objects at a linear resolution limited only 
by the size of the background continuum source. Absorp- 
tion line observations also have a superior sensitivity com- 
pared to the corresponding emission lines and can success- 
fully be used as probes of the composition as well as the 
physical and chemical conditions in the ISM. Such obser- 
vations have for the most part been done at optical wave- 
lengths, but the last few years have seen an increase in the 
number of high redshift absorption line systems observed 
at radio wavelengths. 

The neutral atomic gas component has been observed 
through absorption of the 21cm line of atomic hydrogen. 
About a dozen 21cm HI absorbers at redshifts z > 0.2 
have been detected (Carilli 1995), with typical column 
densities A(HI) = 5 x 10 18 (T sp // H i) cm" 2 (e.g. Briggs 
1988; Carilli 1995), where T sp is the spin temperature in 
K and /hi is the covering factor of atomic gas across the 
extent of the background continuum source. 

Compared to the ionized and atomic parts of the ISM, 
the molecular gas component is characterized by high den- 
sities and low temperatures. Stars are formed directly 
from the molecular gas, and its physical and chemical sta- 
tus gives information about the stellar formation history 
as well as giving the initial conditions for the on-going 
star formation. We have recently detected absorption of 
molecular rotational transitions in four galaxies at red- 
shifts 0.25-0.89 (Wiklind & Combes 1994, 1995, 1996a,b; 
Combes & Wiklind 1996). The background source is in 
all cases a highly obscured radio-loud QSO, located ei- 
ther in the background or in the galaxy where the ab- 
sorption takes place. Absorptions in the millimeter wave 
band provide a velocity resolution more than two orders of 
magnitude higher than in optical spectroscopy. Since the 
amount of molecular gas needed to produce a detectable 
absorption line is less than 1 M Q , the sensitivity of molec- 
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ular absorption line observations in terms of total mass 
is ~ 10 12 times higher than for the corresponding emis- 
sion lines (cf. Wiklind & Combes 1994). This permits us 
to derive physical properties of the molecular component 
of the ISM at high redshifts which would otherwise be 
impossible. 

In this paper we present new molecular absorption 
line data on the absorption system at z=0. 24671 towards 
PKS1413+135. In Sect. 2 we give a general background to 
the radio source PKS1413+135 and the absorbing galaxy. 
The observations and data reduction are presented in 
Sect. 3. Our results are given in Sect. 4: in particular, we 
derive the physical properties of the molecular gas and we 
discuss possible time variations in the absorption lines. In 
Sect. 5 we compare the chemistry of the redshifted molec- 
ular gas with Galactic and nearby extragalactic systems 
and we discuss constraints to variations in physical con- 
stants derived from our absorption line data. Finally in 
Sect. 6 we summarize our results. 

2. PKS1413+135 and the galaxy at z=0.247 

The radio source PKS1413+135 has a flat continuum spec- 
trum extending from a few hundred MHz into the near- 
infrared, with an exponential decrease at ^5 /im (e.g. 
McHardy et al. 1991). The radio loud AGN is invisible at 
optical wavelengths. On the basis of this steep infrared to 
optical continuum, PKS1413+135 was originally classified 
as a red QSO (e.g. Brcgman et al. 1981). The abrupt de- 
crease of the continuum in the near-infrared was modeled 
as a high-energy cut-off in the electron energy distribu- 
tion. Spectroscopic surveys in both the infrared and opti- 
cal bands showed no prominent emission lines (Brcgman et 
al. 1981; McHardy et al. 1994; Perlman et al. 1996). Stocke 
et al. (1992) detected Balmer absorption lines and weak 
[Oil] 3727 A emission at z=0.247. Rapid variability at ra- 
dio and infrared wavelengths, strong and rapid variable in- 
frared polarization and a flat-spectrum radio continuum, 
lead to a classification of PKS1413+135 as a BL Lac ob- 
ject. Subsequent radio imaging using VLBI techniques has 
questioned this classification. The radio source consists of: 
(1) a compact radio core, (2) a jet-like structure on a 
parsec scale, and (3) a counter-jet (Perlman et al. 1996). 
All the features except the compact core have steep spec- 
tral indices, meaning that their contribution to the contin- 
uum flux is neglible at millimeter wavelengths. The pres- 
ence of a counter-jet disagrees with the interpretation of 
PKS1413+135 as a BL Lac object, since these are believed 
to represent radio sources with the jet pointing almost ex- 
actly along the line-of-sight, with the flux boosted by rcl- 
ativistic beaming. The counter-jet would not be visible in 
this scenario. Perlman et al. (1996) model PKS1413+135 
as a typical Compact Symmetric Object (CSO); a young 
radio source with an age < 10 4 yrs. Doppler-boosting of 
the continuum could occur in the core, but not in the par- 
sec scale jet and counter-jet. 



Imaging of the system in optical and infrared bands 
reveals an edge-on galaxy (inclination i ss 87°) with a 
prominent dust lane (McHardy et al. 1991, 1994; Stocke 
et al. 1992). The luminosity profile is exponential with a 
scale length of 1.8", corresponding to 5.8 kpc for H = 
75kms~ 1 Mpc~ 1 and q = 0.5 (these values are used 
throughout this paper). The disk-to-bulge ratio is ~1.1 
(McHardy et al. 1994), but this value is uncertain due to 
unknown extinction corrections. The integrated blue lu- 
minosity (obtained from the V-band which approximately 
corresponds to the B-band in the restframe of the galaxy) 
is 19.6 ±0.15 (McHardy et al. 1992), giving M B = -20.49 
or L B = 2.4 x 10 10 L Q . No sign of the AGN is seen at 
optical wavelengths, but in the infrared the luminosity 
becomes dominated by a point source. This point source 
is situated within 0.1" (325 pc) of the center of the galaxy 
(Stocke et al. 1992). If the radio source is situated behind 
the galaxy, this close alignment would lead to gravitational 
lensing. No signs are found for image multiplicity or dis- 
tortion, in either optical or radio bands (e.g. Stocke et al. 
1992; Perlman 1996), strongly suggesting that the galaxy 
at z=0.247 is host to the radio source. 

The X-ray spectrum shows a deficiency of low-energy 
X-ray photons, indicative of an extinction Ay > 30 mag 
(Stocke et al. 1992; McHardy et al. 1994). This corre- 
sponds to an absorbing screen with N(H)> 10 22 cm~ 2 . 
Carilli et al. (1992) detected redshifted 21cm absorp- 
tion at z=0. 24671, with a column density N(HI)= 1.3 x 
10 19 (T s / '/ 'hi), where T s is the spin-temperature and Jhi 
is the covering factor of atomic gas across the back- 
ground source. Wiklind & Combes (1994) reported molec- 
ular absorption at the redshift of the galaxy, with N(H 2 )= 
4 x 10 20 cm~ 2 . These column densities are much smaller 
than that inferred from the X-ray data, but suggest that 
the galaxy at z=0.247 is of late type. The presence of 
molecular and atomic gas, the exponential luminosity pro- 
file and the broad-band colours are consistent with the 
galaxy being of morphological type Sb-c. 

Hence, PKS1413+135 is a radio source situated in a 
spiral galaxy at a redshift z=0. 24671. The blue luminosity 
of the galaxy is typical for an Sbc. The abrupt cut-off in 
the optical continuum is caused by extinction, and not by 
a cut-off in the electron energy distribution. An excess 
in near-infrared emission should be present due to dust 
heated by radiation from the AGN or by the interaction 
of the relativistic electrons with the ambient interstellar 
medium. No such excess has been observed (cf. Stocke et 
al. 1992). 

3. Observations and data reduction 

3.1. Observations 

The observations presented in this paper have mainly been 
obtained at the IRAM (Institut de Radio Astronomie Mil- 
limetrique) 30-m telescope on Pico Veleta in Spain during 
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Fig. 1. Total averages of the J=0^1 transitions of CO and HCO+, the J=1^2 and J=2^3 transitions of HCN, HCO+ 
and HNC towards PKS1413+135. The velocity resolutions are ~0.2kms _1 for all spectra. All transitions have been obtained 
with the IRAM 30-m telescope, except HCO + (0— >1) which was observed with the NRAO 12-m telescope at Kitt Peak. The 
continuum levels have been normalized to unity and zero velocity corresponds to a heliocentric redshift of 0.24671. 



several observing runs from February 1994 to July 1996. 
Some absorption lines were^lso observed at the NRAO 
12-m telescope at Kitt PeakEl in February 1996. 

At the IRAM 30-m telescope we used 3-, 2- and 1.3- 
mm SIS receivers, tuned to the redshifted frequencies of 
the molecular transitions in question (see Table 1). The ob- 
servations were done with a nutating subreflector, switch- 
ing symmetrically ±90" in azimuth with a frequency of 
0.5 Hz. The continuum levels of the observed sources were 
determined using a continuum backend and increasing the 
subreflector switch frequency to 2 Hz. The image sideband 
rejection of the receivers varied from OdB (double side- 
band) to 20 dB (single sideband), and was taken into ac- 
count to determine the continuum temperatures T c dis- 
played in Table 2. The typical system temperature ranged 
from 200 K and up, depending on the whether conditions. 
Data obtained under bad weather conditions have been 
discarded. For the line observations we used broad filter- 
banks and narrow band autocorrelators, the former giv- 
ing velocity resolutions l-3kms _1 , depending on the ob- 
served frequency, and a bandwidth of 500-1000 km s _1 . 
With the autocorrelators we chosed frequency resolutions 
of 10 kHz for the 3-mm receiver, 20 kHz for the 2-mm 
and 40 kHz for the 1-mm receivers. The corresponding 
velocity resolutions are 30ms _1 , 40ms _1 and 50ms _1 , 
respectively. The bandwidths for the autocorrelators were 
40-50 km s _1 . The pointing of the telescope was checked 
regularly on the observed sources or on nearby continuum 

1 NRAO is a facility of the National Science Foundation, op- 
erated under cooperative agreement by Associated Universities 
Inc. 



sources. Since the weather was good and stable for most 
of the observations, we could monitor the continuum flux 
of the observed source during the integration to directly 
check the pointing of the telescope. Typical pointing cor- 
rections were 5-10" . 

At the NRAO 12-m telescope we used the SIS re- 
ceiver for the lower part of the 3-mm band, tuned to the 
redshifted frequencies of the observed transitions. The re- 
ceiver contains two mixers allowing dual polarization ob- 
servations. For the lowest frequencies, only one of the mix- 
ers could be used due to high receiver temperatures in the 
other. Typical system temperatures were in the range 300- 
500 K. The observations were done with a dual chopping of 
the subreflector; the chopping was done at a rate of 1.25 Hz 
with a beamthrow of 2'. As backends we used the filter- 
banks, giving spectral resolutions of 0.25 and 0.50 MHz. 
The pointing and focus were checked regularly on Venus 
and Saturn. 

The data presented here is given in chopper-wheel cal- 
ibrated antenna temperatures (T£). The NRAO data is 
initially represented in the Tj£ scale. To transform it into 
the Tj£ temperature scale we multiplied the T^j values with 
the forward spillover and scattering efficiency, r/f ss = 0.68. 

3.2. Data reduction 

Each spectrum is the average of several 5-10 minutes inte- 
grations with a calibration in-between. Individual spectra 
were inspected for interference and bad spectral baselines. 
The channel-to-channel noise rms was used as a weighting 
factor when adding individual spectra together. Usually 
no baseline corrections were applied, although a small set 
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3.3. The continuum level 



Table 1. Observed molecules and transitions 



Mol. 


Transition 


^rest 

GHz 


Jb) 

obs 

GHz 


Dot. c > 


Tcl. d > 


CO 


0- 


>1 


115.271204 


92.460319 


Yes 


I 


13 co 


0->l 


110.201370 


88.393748 


No 


1 


HCN 


0- 

1- 

2- 


► 1 

>2 
>3 


88.631847 
177.261111 
265.886180 


71.092593 
142.183115 
213.270271 


No 
Yes 
Yes 


K 


HCO+ 


0- 
1- 
2- 


>1 
►2 
>3 


89.188523 
178.375065 
267.557619 


71.539109 
143.076630 
214.610951 


Yes 
Yes 
Yes 


K 


H 13 CO + 


1- 


>2 


173.506782 


139.171726 


No 


1 


HNC 


1- 
2- 


.2 
>3 


181.324758 
271.981067 


145.442611 
218.159048 


Yes 
Yes 


; 


N 2 H+ 


0- 
1- 
2- 


>1 
,2 

>3 


93.173400 
186.344874 
279.511671 


74.735424 
149.469302 
224.199430 


No 
No 
No 


K 


CN 


0- 


► 1 


113.490980 


91.0323830 


Yes 




CS 


1- 


j. 2 


97.980950 


78.591613 


No 




H 2 CO 




* 2 1,2 


140.839504 


112.968937 


No 




2 


1() -» ll 


118.750 


95.251 


No 





a) The rest— frequency of the observed molecules taken from Lovas 
(1992). 

b) Derived from v obs = u TCBt / (1 + z a ). 

c) Detection: Yes/No. 

d) Telescope: I=IRAM, K=Kitt Peak 12-m. 



of spectra obtained under less favourable weather condi- 
tions were corrected with a 3rd order polynomial. The con- 
tinuum level was derived directly from the spectra, taking 
care to only include data obtained under stable weather 
conditions. At the IRAM telescope we also used a con- 
tinuum backend and a fast chopping technique to derive 
T c . This was done at all the observed frequencies. The 
accuracy of the continuum level is critical for the subse- 
quent interpretation of the data (see below). The spectra 
have not been binned or smoothed except in the cases 
when adding together data obtained with different fre- 
quency resolutions. In these cases spectra are rebinned to 
correspond to the spectra with the lowest spectral reso- 
lution. The spectra presented here have their continuum 
level normalized to unity. 



The observed continuum level T c of the background source 
is an important quantity when deriving parameters from 
the absorption lines. It is, however, difficult to obtain T c 
with an accuracy better than 10% in a single measure- 
ment, due to calibration uncertainties inherent in the at- 
mospheric models used for the chopper-wheel method. 
The variance about the mean of the opacity can be ob- 
tained by differentiating Eq. (4) in Sect. 4.2 and is 



exp(r„) - 1 

°V„ ~ ™ OT c , 



(1) 



where <tt c is the standard deviation of the observed con- 
tinuum temperature and we have neglected the variance 
of T^bs, the measured depth of the absorption line, which 
generally is much smaller than that of the continuum tem- 
perature. For t « 1, <t Ti ,/t w <7t c /T c , which is usu- 
ally small for finite values of r and T c . For r >> 1, 
<j Tv « exp(T)CTT c /r c and we loose all information about 
the true opacity. In the critical opacity region (r il 1), the 
error in our estimate of the opacity depends nonlinearly 
on the continuum level and the opacity itself. 

More specifically, the important parameter when de- 
riving the opacity of an absorption line seen towards an 
unresolved background source, is the line-to-continuum 
ratio. The values derived are therefore independent of the 
true value of the continuum level as long as the only vari- 
ations in T c originate in instantaneous pointing uncer- 
tainties. This is the situation under good weather con- 
ditions, where the average calibration uncertainty is the 
only parameter affecting the derived opacity. When the 
atmospheric conditions change between two calibrations 
(typically 5-10 minutes), we loose information about the 
continuum level and the line-to-continuum ratio ceases to 
be a good measure of the opacity. When this has been the 
case, we chose a period of acceptable atmospheric condi- 
tions and adopted the derived T c from this period to scale 
the rest of the data. However, in this procedure we be- 
come sensitive to pointing errors. Fortunately, their vari- 
ance was always lower than 20% of the beam-width. 

The main uncertainty in the derived column densities, 
excitation temperatures and, to a lesser extent, abundance 
ratios, originates in the uncertainty of the continuum level. 
The nonlinear dependence makes it difficult to quantify 
the formal uncertainty. Formal errors presented here are 
therefore based on the assumption that T c is correct. The 
relative errors in T c are then to be multiplied by 1.3, 1.7 
and 3.2 to obtain the relative errors induced on r, when 
t v = 0.5, 1 and 2 respectively. 

4. Results 

4-1. The observed absorption lines 

Towards the radio source PKS1413+135 we have observed 
18 different molecular rotational transitions of a total of 
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11 molecules. Of these we clearly detect 9 transitions of 5 
different molecules; the J=0— >1 lines of CO, HCO + and 
CNl, the J=l->-2 and J=2->3 lines of HCN, HCO+ and 
HNC. The observed transitions with their rest-frequencies 
and corresponding rcdshifted frequencies are presented in 
Table 1. In Table 2 we list the observed integrated opti- 
cal depths and the derived excitation temperatures and 
column densities. Notice that some of the column densi- 
ties in Table 2 are modified after a more thorough analysis 
(Sect. 4.2.3). The final column densities are presented sep- 
arately in Table 3. 

The detected transitions are shown in Figs. 1 and 2. 
The continuum levels have been normalized to unity. All 
the spectra in Fig. 1 except HCO + (0^1) are averages of 
the observing sessions in 1994 (cf. Sect. 3). The 1995 and 
1996 data have not been used since there might be tem- 
poral changes in the absorption profiles (see Sect. 4.5). 
In order to show the weak CN(0— >1) line, the scaling 
in Fig. 2 is different from that of Fig. 1. The CN hyper- 
fine transitions which fall within the spectrometer are 
the J= 3/2 - 1/2 with F= 1/2 - 3/2, F= 3/2 - 3/2, 
F= 1/2 - 1/2, F= 5/2 - 3/2 and F= 3/2 - 1/2. Only the 
F= 5/2 — 3/2 line is formally detected, possibly blended 
with the F= 3/2 — 1/2 line. The height of the markings 
corresponds to the observed intensities in the Orion hot 
core (cf. Lovas 1992). 

Some of our non-detected transitions are shown 
in Fig. 3. Since the isotopic lines of 13 CO(0^1) and 
H 13 CO + (1^2) remain undetected at low noise levels, the 
corresponding 12 CO and H 12 CO + lines are not saturated, 
despite absorption depths close to unity. This means that 
the derived column densities are not affected by satu- 
ration effects. The peak opacities for the CO (0— >1) and 
HCO+(l->2) lines are 2-3 (see Table 2). 

The column density of CO is estimated to be ap- 
proximately 2 x 10 16 cm~ 2 . With a Galactic CO to H 2 
abundance ratio of 5 x 1CP 5 (e.g. van Dishoeck & Black 
1987), the column density of molecular hydrogen towards 
PKS1413+135 at z=0.24671 is 



4-2. Derivation of physical parameters 

In order to derive physical parameters from molecular ab- 
sorption line data one must consider three uncertainties: 
(i) the covering factor of absorbing gas across the finite 
extent of the background continuum source, (ii) the exci- 
tation temperature of the molecular gas and (Hi) if LTE 
conditions apply. The first uncertainty affects the deter- 
mination of the velocity integrated optical depth while the 
second and the third also affects the derivation of column 
densities. 
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Fig. 2. Normalized spectra of CN N=0— >1 lines. From left 
to right we have indicated the J= 3/2 - 1/2, F= 1/2 - 3/2, 
F= 3/2 - 3/2, F= 1/2 - 1/2, F= 5/2 - 3/2 and F= 3/2 - 1/2 
hyperfine transitions. Only the F= 5/2 — 3/2 line is formally 
detected, possibly blended with the F= 3/2 — 1/2 line. The 
height of the markings corresponds to the observed intensities 
in the Orion hot core (cf. Lovas 1992). The velocity resolution 
is 3.3kms _1 and zero velocity corresponds to a heliocentric 
redshift of 0.2467. 



4.2.1. Opacity 

In a manner similar to Wiklind & Combes (1994) we ex- 
press the observed continuum temperature, T c , away from 
an absorption line as 



T c = f s J(T b ), 



(2) 



2 For CN we clearly detect one hyperfine transition: N=l-0, 
J=3/2-l/2, F=5/2-3/2. 



where / s is the beam filling factor of the region emit- 
ting continuum radiation, Tb is the brightness tempera- 
ture of the background source and J(T) — (hv/k)/[\ — 
exp(— hf/kT)]. The spatial extent of the region emit- 
ting continuum radiation at millimeter wavelengths is un- 
known. It is generally accepted that as the wavelength 
decreases, the emission originates from regions closer to 
the central AGN. At centimeter wavelengths, the core 
of PKS1413+135 is unresolved at milliarcsecond (mas) 
scales (Perlman et al. 1996). The core is likely to be 
smaller at millimeter wavelengths. The BL Lac 3C446 
has been observed with mm-VLBI and has an extent 
< 30 /larcseconds (Lerner et al. 1993). Since the angular 
size of the telescope beam is typically 25", the bright- 
ness temperature of the background source Tb is at least 
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Table 2. Observed properties for PKS1413+135 



Molecule 


Transition 


c 
mK 


f 

kms- 1 


V 


c) 
"t 

km," 1 


s v d 

km," 1 


Date(T) 




X 

K 




N 


2 




CO 


0— >1 
0— >1 
0— >1 
0— >1 
0— >1 


1000 

240 
495 
338 
230 


3.64 ± 0.06 

3.07 ± 0.11 
3.62 ± 0.12 
3.78 ± 0.08 

3.08 ± 0.36 


2.8 
2.5 
>2.0 
>2.4 
>1.4 


0.063 
0.077 
0.138 
0.097 
0.290 


0.06 
0.13 
0.06 
0.06 
0.13 


Total 

Apr95(I) 

May94(I) 

Mar94(I) 

Feb94(I) 








„ „ + 0.8 
2d -0.7 


X 


10 16 


13 co 


0— >1 
0— >1 
0— >1 


1000 
240 
512 


<0.14 
<0.23 
<0.19 


<0.04 
<0.06 
<0.07 


0.043 
0.019 
0.022 


0.10 
1.06 
0.72 


Total 

Apr95(I) 

May94(I) 






< 




X 


10 14 


HCN 


0— >1 


75 


<0.92 


0.06 


0.063 


2.08 


Feb96(K) 
















1— >2 
1^2 
1— »2 
1— »2 


1000 
170 

320 
243 


0.99 ± 0.07 
1.52 ± 0.14 
1.36 ± 0.11 
0.46 ± 0.10 


0.5 
0.6 
0.6 
0.4 


0.042 
0.077 
0.079 
0.079 


0.21 
0.21 
0.16 
0.24 


Total 

Apr95(I) 
May94(I) 
Mar94(I) 


11 


. 8± 3.1 






X 


10 12 




3-2 
2-»3 
2 -.3 
2-»3 


1000 
120 
173 
146 


0.93 ± 0.10 
0.95 ± 0.53 
1.53 ± 0.19 
0.80 ± 0.10 


0.6 
0.4 
0.8 
0.6 


0.063 
0.233 
0.135 
0.085 


0.22 
0.33 
0.16 
0.16 


Total 

Apr95(I) 
May94(I) 
Mar94(I) 


11 


=+3.1 
8 -2.1 




7 3+ 11 
7d -0.3 


X 


10 12 


HCO + 


0-»l 


46 


5.20 ± 0.50 


1.5 


0.218 


0.41 


Feb96(K) 


5 


4+ - 7 
-0.5 












1— >2 
1— >2 
1— >2 
1— >2 
1— >2 


1000 
170 
340 
242 
128 


6.87 ± 0.06 

5.57 ± 0.19 
7.22 ± 0.11 
6.20 ± 0.15 
4.51 ± 0.32 


>2.6 
>2.0 
>2.7 
>2.3 
>1.4 


0.079 
0.131 
0.070 
0.108 
0.268 


0.04 
0.21 
0.12 
0.16 
0.12 


Total 

Apr95(I) 
May94(I) 
Mar94(I) 
Feb94(I) 


8 


7 +0.2 
-0.2 




29 +o.i 

-0.2 


X 


10 13 




2^3 
2^3 
2— >3 


1000 
120 
172 


4.75 ± 0.10 
5.98 ± 0.19 
4.01 ± 0.18 

A 70 -U n 

'i. 1 z ztz U . o I 


2.0 
2.3 
1.3 


0.122 
0.085 
0.172 


0.06 
0.32 
0.11 


Total 

Apr95(I) 
May94(I) 

IVTm-O/1 ( ~i \ 

IVltlir :J 4t { I ) 


8 


,7 io.2 


o g+0.05 
29 -0.05 


X 


10 13 


H 13 CO + 


1— >2 


329 


<0.18 


<0.08 


0.027 


0.42 


May94(I) 


8 


7 +0.2 
-0.2 


< 


7 9 + - 1 
'■' -0.2 


X 


10 11 


HNC 


1— >2 
1— »2 
1— »2 
1— »2 


1000 
170 
310 
238 


0.53 ± 0.05 
0.68 ± 0.11 
0.69 ± 0.09 
0.45 ± 0.07 


0.3 
0.3 
0.3 
0.2 


0.040 
0.043 
0.078 
0.054 


0.21 
0.41 
0.12 
0.16 


Total 
Apr95(I) 
May94(I) 
Mar94(I) 


15 


. + 9.5 
-5.0 




5 2 + 5 - 7 


X 


10 12 




2-»3 
2^3 
2^3 


1000 

120 
168 


0.60 ± 0.11 
0.93 ± 0.29 
0.78 ± 0.18 


0.3 
0.4 


0.085 
0.127 
0.114 


0.21 
0.32 
0.21 


Total 

Apr95(I) 
May94(I) 


15 


■^ 9 5: 6 o 




5 2 + 3 - 2 


X 


10 12 


CS 


1— >2 


80 


< 0.62 


<0.21 


0.048 


1.88 


Feb96(K) 






< 


2 . +0.7 


x 


10 13 


CN b ) 


0-»l 
0— >1 
0^ 1 
0— >1 


1000 
240 
341 
150 


1.10 ± 0.23 

< 0.19 

< 0. 15 
1.10 ± 0.23 


0.1 
<0.03 
<0.01 

0.1 


0.012 
0.031 
0.014 
0.012 


3.29 
1.03 
3.29 
3.29 


Total 

Apr95(I) 
Mar94(I) 
Jul96(I) 








2.o+; : 2 


X 


10 13 


N 2 H + 


0-»l 


075 


< 1.07 


<0.10 


0.103 


1.00 


Feb96(K) 






< 


q q+3.7 


X 


10 12 




1— »2 


170 


< 0.48 


<0.03 


0.028 


2.00 


Apr95(I) 




io+J 2 


< 


2 1 + 0.5 
-0.4 


X 


10 12 




2^3 


120 


< 1-1 


<0.08 


0.077 


1.34 


Apr95(I) 




ro+_l 


< 




X 


10 12 


H 2 CO 


1 1,1 — * 2 1,2 


296 


< 0.14 


<0.03 


0.031 


0.55 


Mar94(I) 




u,+l 


< 




X 


10 13 


o 2 


lO — 1 ! 


484 


< 0.14 


<0.03 


0.031 


0.55 


May94(I) 




io+l 


< 


2 . +0.7 


X 


10 16 



a) The continuum level in mK. Total averages have T c = 1000 mK. 

b) The maximum optical depth. A lower limit is set when the absorption reaches within the spectral noise rms of the zero level. 

c) Channel— to— channel noise rms of the opacity. 

d) Channel separation used for the derivation of opacities. 

e) Derived or assumed excitation temperature in K. 

f) Derived for an upper limit to the N= 1-0, .1= 3/2 - 1/2, F= 5/2 - 3/2 lino. 
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opacity t v of the gas. AT£ = T c — |T a b s |, where T a b s is the 
depth of the absorption line measured from the continuum 
level. When r v is very large the absorption saturates, i.e. 
T* hs — ► T c and AT£ — > 0. In this equation, we cannot 
derive independently the covering factor of the absorbing 
clouds, f c , and their optical depth. Combining Eqs (1) and 
(2) we can express the 'true' optical depth as: 



1 



t v = - ln[l - T (l - cxp(-T„ ob J)] 

Jc 



where 



(4) 



T ^ob S ~ ^ AT* ^ ' 



T„ oba is the optical depth averaged over the size of the 
background continuum source and over the excitation con- 
ditions of the molecular gas along the line of sight. This 
value is given for each observed absorption line in Table 2. 
By assuming that the opacity of the absorbing gas is high, 
Eq. (3) gives a minimum possible value for the covering 
factor f c 1 — exp(— T„ obs ). Should the opacity be low, the 
filling factor must be even larger. The velocity integrated 
opacity is obtained by integrating t v in Eq.(4) over the 
full extent of the absorption line. 

In PKS1413+135 the optical depth of the CO(0->1) 
and HCO+(1^2) lines is ~3. This means that the ab- 
sorption almost reaches the zero level, but it is not heavily 
saturated. The HCN, HNC and the J=2^3 line of HCO+ 
all have smaller optical depths. The non-detections of less 
abundant molecules such as CS, N2H+ and the isotopic 
variants 13 CO and H 13 CO + are consistent with the de- 
tected lines not being heavily saturated. The large depth 
of the CO(0->1) and HCO+(1^2) lines then implies that 
the covering factor / c must be very close to unity. In the 
following we will set f c = 1 for PKS1413+135. 



Fig. 3. Normalized spectra of the non-detected 13 CO(1-0), 
H 13 CO+(1^2) and O 2 (l -> li) lines. The velocity resolution 
is 0.65, 0.66 and 0.64kms~ , respectively. The velocity scale is 
heliocentric. Notice the different scale for the continuum level 
compared with Fig. 1. 



10 9 x T c . This means that we do not have to consider the 
local excitation temperature of the molecular gas when 
deriving the opacity, since this is typically of the order 
10 K. The excitation does enter, however, when deriving 
column densitites. 

The observed antenna temperature of an absorption 
line can be expressed as 



ATXH = f s J(T h ){l - f c [l - exp(-r v ))}, 



(3) 



where f c is the fraction of the continuum source area cov- 
ered by molecular gas at the velocity corresponding to the 
frequency v. We here assume that f c is a characteristic 
covering factor for the molecular gas, independent of the 



4.2.2. Excitation temperature and column density 

The excitation temperature T x relates the relative popu- 
lation of two rotational levels of a molecule as: 



— = — exp (-hv 2 i/kT^) , 



(5) 



where gi is the statistical weight for level i and hv 2 \ is 
the energy difference between two rotational levels. In 
order to derive T x we must link the fractional popula- 
tion in level i to the total abundance (here abundance is 
used in synonym with column density). This is done by 
invoking a weak LTE-approximation and assuming that 
T x = T rot , where T rot is a temperature which governs the 
fractional population of all rotational levels in a given 
molecule. The LTE approximation is weak in the sense 
that it does not imply that T rot equals the kinetic temper- 
ature and it allows for different molecules to have different 
T r ot- With this approximation we can use the partition 
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Fig. 4. a) The function /(T x ) for the 4 lowest rotational tran- 
sitions of HCO + as a function of the excitation temperature 
T x (cf. Eq. 6) . b) The ratio of integrated optical depth for the 
J=1^2 line of HCO + for two different excitation temperatures, 
Ti and Ti. The number at each line in the plot corresponds to 
the factor a, Tb = aTi. The results shows that for a given tran- 
sition and equal column densities, the coldest gas component 
contributes the largest part of the observed opacity. 
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Fig. 5. a) A schematic illustration of a two-component molec- 
ular medium seen in absorption. The two components are char- 
acterized by excitation temperatures Ti and Tb, with Tb > Ti 
and represents a fraction 7 and 1 — 7 of the total column den- 
sity N (see the text for more details); b) Same as a) but with 
the warmer (denser) gas distributed in a clumpy medium. The 
line of sight is shown 'extended' to emphasize the dumpiness 
present within the gas causing the absorption profile. 



function Q(T X ) — J2y =0 gjexp(—Ej/kT K ) to express the 
total column density N to t as 



Ntot = 



8tt 



c 3 g.jAj J+1 

Q(T x )exp(Ej/kT x ) 
1 - exp(-/u//fcT x ) 



7 x) f r v dV , 



(6) 



where / r v dV is the observed optical depth integrated over 
the line for a given transition, gj = 2 J + 3 for a transi- 
tion J — > J + 1 and _Ej is the energy of the rotational 
level J. By taking the ratio of two observed transitions 
from the same molecule, the excitation temperature can 
be derived. If the lower transition is saturated, T x will be 
an upper limit. Non-saturated lines will directly give T x , 
limited only by the S/N in the spectra. The derivation of 
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Fig. 6. The 
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surface, of the 
PKS1413+135 
J=1^2 and J=2->3 transitions of HCO+. 
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the excitation temperature depends only on the line-to- 
continuum ratio of the two transitions in question and is 
therefore independent of a correctly measured continuum 
level, as long is this is not caused by changing weather con- 
ditions inbetween two calibrations. The T x derived in this 
way is an average excitation temperature along the line of 
sight and averaged over the finite extent of the continuum 
source. The strong frequency dependence of the column 
density in Eq,(6) is only apparent since the Einstein coef- 
ficient is proportional to v 3 . 

4-3. Multicomponent analysis 

Towards PKS1413+135 we have observed three transitions 
of HCO+ which allows us to derive the excitation tem- 
perature using the ratio of the observed integrated opaci- 
ties for the J=0^1/J=1^2 and J=1^2/J=2^3 lines. If 
weak-LTE conditions apply, as defined above, these ratios 
should give the same excitation temperature. This, how- 
ever, is not the case. The first ratio gives T x = 5.4+g'g K, 
the second 8.7t ' 2 K. Is this an indicaton that the assump- 
tion of LTE conditions are at fault? 

The above derivation of the excitation temperature 
and column density corresponds to values averaged over 
the angular extent of the background source as well as 
along the line of sight. Since these extents are both finite, 
we may well have several gas components characterized 
by different excitation temperatures contributing to the 
observed absorption lines. 

The observed property of a molecular absorption line is 
the velocity integrated optical depth J T„dV, which is pro- 
portional to iV to t//(T x ) (e.g. Eq. 6). As shown in Fig. 4a, 
the function f(T K ) is highly nonlinear. The ratio of the in- 
tegrated opacities for a given molecule and transition will 
therefore depend strongly on the excitation temperature. 
The ratio of opacities of two gas components characterized 
by excitation temperatures T\ < T 2 is shown in Fig. 4b for 
HCO+(1^2), for T 2 /T 1 ratios of 2-10. It is evident from 
the figure that the 'cool' gas component contributes sig- 
nificantly more to the opacity than the 'warm' gas. Note 
that the temperature refers to the excitation conditions, 
an excitationally cold gas can still be kinetically warm if 
the density is low enough. It is thus possible to combine 
a relatively large column of warm (i.e. dense) molecular 
gas with smaller column of excitationally cold (i.e. diffuse) 
gas to produce the observed opacity. In order to see how 
such a two-component medium can combine and still be 
consistent with the observed data, we express the opacity 
as: 



Fig. 7. The same fraction 7 as in Fig. 6, but now including 
solutions also for the J=0— »1/J=1— 1-2 ratio. The grey-scaled 
area in the lower part of the figure represents the region where 
the two solutions agree, suggesting 7 = 0.2, T CO oi = 4.2 K and 
Twarm ~ 25 - 30 K. 



J T v dV = 



8tt 



7 



1-7 

m) ' f(T 2 ) 



+ 



(7) 



where 7 is the fraction of the total column density char- 
acterized by a temperature T\ and 1 — 7 is the fraction 
associated with gas of a temperature T 2 . In the following 
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we will assume that T\ < T 2 and that the covering fac- 
tor of the two components are equal. A schematic view 
of this two-component model is shown in Fig. 5a. In the 
next section we will study the case when the assumption 
of equal filling factors is relaxed. Taking the ratio of two 
observed opacities the column density is eliminated and 
we can solve for 7 as a function of 7\ and T 2 . One result 
obtained for the J=l — >2 and J=2^3 lines of HCO + is 
shown in Fig. 6. The ratio of the observed opacities give 
an average excitation temperature of 8.7 K (see Table 2), 
which is the limit for 7 = 1 in Fig. 6. It is, however, pos- 
sible to have a fraction 7 <1 of the total column density 
in the form of significantly colder gas, with the remaining 
1 — 7 in the form of a warmer gas component. The a;— axis 
in Fig. 6 shows the ratio of the two temperatures. Hence, 
one solution is to have 25% of the column density in the 
form of a gas component with T x «5K and 75% in the 
form of a warm gas component with T x «20K. 

By including all three transitions observed for HCO + , 
we get two sets of solutions for the fraction 7. These are 
shown in Fig. 7. The grey-scaled area indicates where the 
two solutions coincide. The results suggests the presence 
of an excitationally cold component with T x « 4 K which 
constitutes a fraction ~20% of the total column density. 
The remaining 80% is made up of a gas with a T x w 20 — 
40 K. Hence, the non-concordant values of the excitation 
temperatures derived from three transitions of HCO + do 
not necessarily imply that the assumption of weak-LTE is 
violated, but are consistent with the presence of two gas 
components with different excitation temperatures. 

There might be more than two gas components along 
the line of sight. By making the ansatz that we have N 
components, each charaterized by an excitation tempera- 
ture Ti, with Ti < T 2 < ... < T N , the fraction 7 of the 
coldest gas component will be somewhat smaller than de- 
rived assuming a 2-component model. However, the strong 
nonlinear behaviour of the opacity as a function of T x 
(Fig. 4a, b), ensures that a third, fourth or higher gas com- 
ponent will give a small or neglible contribution to the ob- 
served opacity. Thus, a 2-component model should give a 
reasonable approximation to the possible gas components 
existing along the line of sight and over the finite extent 
of the background source. 

Adopting values for T\ and T 2 , giving a fraction 7, the 
total column density N can be derived. This value is al- 
ways larger than for the 'average' T x , with a maximum 
for the smallest allowed fraction 7 of cold gas. This is a 
natural consequence of the nonlinear dependence of the 
opacity to the excitation temperature as shown in Fig. 4a. 
In the case of HCO + , the allowed solutions give total col- 
umn densities in the range 3.7 — 11.2 x 10 13 cm~ 2 , with 
corresponding excitation temperatures 3.4 & 13.3 K and 
4.5 & 44.5 K for the cold and warm components. The frac- 
tion 7 is 0.19 and 0.12 for these two cases. The column 
density derived using the average T x as given in Table 2 
is 2.9 x 10 13 cm" 2 . The maximum allowed fraction of cold 



gas is 7 = 0.22, with excitation temperatures 3.9 & 18.0 K 
and a total column density of HCO+ of 4.5 x 10 13 cm~ 2 . 

4-4- A clumpy medium 

Since the excitation of molecular gas is generally governed 
by collisions with H 2 , a high-X^ gas is characterized by a 
high density.The division in cool and warm gas made in 
the previous section is therefore equivalent to a division 
into diffuse and dense gas. A more realistic multi-phase 
medium combining the dense and diffuse phases is that 
observed in the Milky Way ISM, where a clumpy struc- 
ture of small and dense clouds with a small filling factor 
is embedded in a pervasive diffuse molecular medium. A 
schematic view of this two-component model is shown in 
Fig. 5b. 

If we call / c the area filling factor of each clump, for 
which t„ >> 1, their contribution to the absorption depth 
is simply f c (assuming that the brightness temperature of 
the background source is homogeneous over its extent). If 
we assume that there are in total n c clouds in the beam, 
and that their line-width is 5v each, corresponding to a 
velocity filling factor /„ = Sv/AV (where AV is the to- 
tal spectrum line- width observed), and neglecting cloud 
overlapping for the sake of simplicity, their combined con- 
tribution to the opacity, at at given v, will be 



n-cfvfc = 1 - CXp(-T bs) , 



(8) 



where r b s is the observed opacity. The total H 2 column 
density in the beam is then 



^H 2 



-nB_ 2 r c - 



1 - exp(-T obs ) 

fcfv 



cm 



(9) 



where tih 2 is the density of molecular hydrogen in each 
clump and r c the characteristic radius of the clumps. 
In order to see how the presence of small dense struc- 
tures influence the derivation of column density we can 
use the simple parametrization presented above, together 
with the characteristic values derived from observations 
of the Galactic molecular ISM. Assuming a r b s = 0.5 
for the CO(0— +1) line a velocity width Sv — 0.1 kms -1 
and a covering factor f c = 0.001 per clump, the num- 
ber of small clumps in the beam will be n c w 4 x 10 3 . 
If we further assume a radius of 100 AU for the clumps 
(approximately corresponding to 0.1% of the area of a 
background source with a diameter of 6000 AU) and a 
density nn 2 = 10 4 cm~ 3 , Eq. (10) gives a column density 
N{H 2 ) ~ 8 x 10 22 cm~ 2 . The corresponding CO column 
density, derived solely from the observed opacity and as- 
suming f c « 1, would in this case be 3 x 10 15 cm~ 2 , which 
translates to a H 2 column density of 6 x 10 19 cm -2 . This 
simple excercise shows that the column density can be un- 
derestimated by three orders of magnitude if the molecular 
ISM is made up of small dense clumps rather than a dif- 
fuse medium. 
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Fig. 8. a) Normalized spectra of the J=0^1 transition of CO obtained in March 1994, April 1995 and July 1996 and the 
J= 1— >2 transition of HCO + obtained in March 1994 and April 1995. The velocity resolution is 0.25 km s -1 for the CO spectra 
and 0.20 kms -1 for the HCO + spectra, b) Normalized spectra of the J=0^1 transition of CO and the J=l— >2 transition of 
HCO + both obtained in May 1994 with a velocity resolution of 0.095 km s _1 and 0.082 km s -1 , respectively. The velocity scale 
is heliocentric and the dashed line corresponds to zero velocity for a redshift z = 0.24671. 



4-5. Time variability 

In Fig. 8 we show the CO(0->1) and HCO+(1^2) absorp- 
tion lines for different epochs, spanning 2 years for CO 
and 1 year for HCO + . In Fig. 8b we show the CO and 
HCO + line at higher velocity resolution. While there are 
no apparent changes within the noise in the HCO + spectra 
between March 1994 and April 1995, the CO line appears 
to have changed. A closer inspection of the CO absorption 
reveals that it is the deepest absorption component that 
has changed while the second deepest component remains 
more or less constant. In order to quantify possible vari- 
ations we fit 3 Gaussian components to the CO spectra 
from all observed epochs. The components are illustrated 
in Fig. 9 and the Gaussian parameters are listed in Table 4, 
together with the integrated opacities derived from the 



Gaussian fits. The data from February 1994 are of poor 
quality, as seen from the high x 2 value of the Gaussian 
fit, and will not be used in the analysis. The integrated 
opacities show that while the second deepest absorption 
component (no. 2 in Table 4) does not show any significant 
variation, the main component (no. 1) changes by a factor 
of > 2 in velocity integrated opacity. The third component 
is much weaker and less clearly defined. The Gaussian fit 
for this component is therefore more questionable than for 
components 1 & 2, which show consistent center velocities 
and half widths for the different epochs. 

How certain are these results? Changes in the inte- 
grated opacity could be mimicked by an incorrectly de- 
rived continuum level, if this is caused by calibration er- 
rors. However, pure pointing uncertainties would not af- 
fect the ratio of the components. It is difficult to quan- 
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Table 4. Gaussfits to the CO(O-fl) line. 



Date 



Comp. 



km s 



km s 



1/2 
1 



[I 



r u dV\ 



km s 



[I 



T„dVY R g) 

J tot 
km s~ 1 



Feb94 


0.46 


1 


0.91 


-0.7 


0.9 






2 


0.38 


+0.5 


1.9 






3 


0.07 


+ 1.1 


7.4 


Mar94 


0.68 


1 


0.97 


-0.5 


1.3 






2 


0.56 


+0.7 


0.6 






3 


0.25 


+ 1.8 


1.1 


May94 


1.00 


1 


1.10 


-0.6 


1.2 






2 


0.47 


+0.6 


0.8 






3 


0.15 


+ 1.6 


1.9 


Apr95 


0.48 


1 


0.83 


-0.6 


1.1 






2 


0.62 


+0.6 


0.6 






3 


0.20 


+ 1.7 


1.9 


Jul96 


0.31 


1 


0.77 


-0.6 


1.5 






2 


0.56 


+0.8 


0.8 






3 


0.31 


+ 1.8 


0.6 



4.87 



0.44 



0.83 



0.31 



0.60 



2.6 
1.5 
0.9 

2.8 

0.48 
0.31 

>3.3 
0.47 
0.32 

1.6 

0.49 
0.42 

1.9 

0.61 
0.63 



3.2 ± 0.2 



3.80 + 0.05 5.6 



4.07 + 0.07 > 7.0 



+0.8 



+0 

0.! 



3.00 ± 0.06 3.3+ ( J- 3 



2.98 + 0.06 3.1 



+0.2 



a) Continuum flux normalized to the flux at May 1994. 

b) The depth of the absorption line from the normalized continuum level. 

c) The center velocity of the Gauss component. 

d) Full width at half intensity. 

e) Integrated opacity of individual components derived from the Gaussian fits. 

f) Total integrated opacities from the observed profiles. 

g) Ratio of integrated opacities of component 1 and 2. Errors correspond to ±10% variation of the continuum level. 



Table 3. Column densities 



Molecule 




^cold, diffuse 


^warm.dcnsc 


N 






K 


K 


cm~ 2 


CO 


0.22 


3.9 


18.0 


2.0 x 10 16 


HCN 


0.05 


3.4 


27.0 


2.0 x 10 13 


HNC 


0.05 


3.4 


27.0 


1.0 x 10 13 


HCO+ 


0.22 


3.9 


18.0 


4.5 x 10 13 



a) The fraction of the column density in the cold gas 
component. 



tify the errors associated with multicomponent Gaussian 
fits to absorption spectra. The errors associated with the 
total velocity integrated opacity are small (see Table 4). 
Hence, as long as the Gaussian fit is good, the error in the 
integrated opacity for the separate components is small 
(^5%). The dominating uncertainty is associated with 



errors in the continuum level caused by an unstable at- 
mosphere. This can be minimized by taking the ratio of 
integrated opacities for the first and second component 
(last column in Table 4), but due to the non-linear rela- 
tion between the depth of the absorption profile and the 
opacity a small uncertainty remains. The continuum level 
is derived using <;10 individual spectra plus a number of 
special continuum observations. Each of these has a typi- 
cal uncertainty of ±15%. The final uncertainty associated 
with the continuum level from the average of these individ- 
ual measurement is approximately ±5%. In order to quan- 
tify the uncertainty of the integrated opacities introduced 
by errors in the continuum level, we have done the analysis 
(Gaussian fits and derivation of opacities) with a contin- 
uum level varying ±10% around the observed value. This 
is a conservative estimate and nevertheless shows that the 
ratio of the integrated opacity of the first and second ab- 
sorption components has varied from <;7 in May 1994 to 
-3 in 1995 and 1996. The ratio of R between May 1994 
and July 1996, including the assumed errors of ±10% in 
the continuum level, is 2.3 ± 0.3. 

The fact that we do not see any clear changes in the 
HCO + line is not contradictory. The opacity of the HCO + 
absorption is higher than that of CO, making it less sen- 
sitive to changes in the continuum level. Moreover, the 
abundance ratio HCO + /CO is higher in low density envi- 
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Fig. 9. Three component Gaussian fitted to the CO(0— »1) ab- 
sorption for the spectrum obtained in March 1994. The residual 
is shown above the spectrum. The parameters of the Gauss fit 
is given in Table 3. The numbers correspond to the components 
as listed in Table 3. 

ronments than at high densities (cf. Lucas & Liszt 1996), 
with the HCO + molecule existing in low density regions 
where CO self-shielding is not taking place. 

5. Discussion 

5.1. Small scale structure in the molecular gas 

In the Milky Way, several lines of evidence suggest that 
the molecular ISM is clumpy on all scales from 20 AU to 
100 pc and characterized by a fractal structure (cf. Fal- 
garone et al. 1992). Time variations in Galactic H 2 CO ab- 
sorption towards extragalactic point sources (Marscher et 
al 1993) even suggest structures in the molecular ISM on 
scales of 10 AU. Through numerical simulations Marscher 
& Stone (1994) were able to derive constraints on the frac- 
tal structure of molecular clouds from the time variability 
detections. The mean number of small clumps along the 
line of sight should be larger than previously thought, i.e. 
the size spectrum of clumps should be a steeper power- 
law, constraining the fractal dimension. Detailed multi- 
level studies of CO and CS near the edge of a nearby 
molecular cloud complex have revealed that the molecu- 
lar ISM is confined to small structures (down to 35 AU) 
which are dense (nn 2 ~ 10 4 cm -2 ), cold (Tk ^10—15 K) 
(Falgarone & Phillips 1996) and characterized by a frac- 
tal structure (Falgarone et al. 1992). In fact, Falgarone 



& Phillips (1996) find no evidence for a pervasive diffuse 
molecular ISM, although this could be a peculiarity for 
the specific region studied. Although emission studies are 
not sensitive to the pervasive diffuse medium, the absorp- 
tion studies of Lucas & Liszt (1996) and Liszt & Lucas 
(1996) show a higher rate of molecular absorption in the 
Galaxy than was originally expected, and precisely detect 
the low-opacity, large-covering factor medium. 

Small scale structures are also present in the atomic 
ISM. Multi-epoch observations of 21 cm absorption 
against high velocity pulsars have detected opacity varia- 
tions of the ISM on a range of scales from 5 AU to 100 AU 
(Frail et al. 1994). VLBI observations at 21cm in front of 
3C radio sources have shown clumps on a scale of 25 AU 
(Diamond et al. 1989, Davis et al. 1996). In the case of 
21 cm absorption against pulsars, where small-scale opac- 
ity structures are detected towards all line of sights, the 
mean opacities are between 0.1 and 2.5, corresponding to 
N(HI) as low as 10 19 - 10 20 cm~ 2 (Frail et al. 1994). The 
opacity variations dectected are quite high (half of them 
have St > 0.1, Frail et al. 1994), so they cannot be ac- 
counted for by mild density fluctuations. Clumps with 
density larger than 10 5 — 10 6 cm~ 3 are implied (Moore 
& Marscher 1995). These clumps could represent 10-20% 
of the total column density. Small dense clumps in the 
atomic gas has also been inferred in high redshift radio 
galaxies (van Oijk et al. 1997), where absorption within 
the Ly-a emission regions suggest the presence of neutral 
atomic clumps of size ^0.03 pc. 

Is the molecular gas along the line-of-sight towards 
PKS1413+135 similar to that of our Milky Way concern- 
ing small scale structure? Combining our molecular ab- 
sorption line data with X-ray data gives three arguments 
which suggest the presence of a molecular ISM consisting 
of small dense clumps embedded in a more diffuse medium: 

— The excitation temperatures for three transitions of 
HCO + arc consistent with a multicomponent molecu- 
lar gas. If two or more components coexist along the 
line-of-sight, the assumption of a single gas compo- 
nent underestimates the true column density. 

— The total column density averaged over the extent 
of the background source and along the line-of-sight 
(atomic and molecular) is much lower than that de- 
rived from the deficiency of low energy X-ray photons. 
Since the covering factor of gas appears to be high, this 
implies that we indeed underestimate the total column 
of gas. 

— Possible time variations in the depth of the molecu- 
lar absorption lines indirectly suggests the presence of 
small scale structure. 

The first item has been presented in Sect. 4.3 and 4.4. 
Below we discuss the second and third item in more detail. 
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Fig. 10. Column density of HCO + plotted versus column densitites of HCN, HNC and CS. The open circles represent Galactic 
diffuse clouds (Lucas & Liszt 1994, 1996), filled circles represent data from Cen A (Wiklind & Combes 1997) and open tri- 
angles represent absorption data towards SgrB2 (Greaves & Nyman 1996). The filled star represents our absorption data for 
PKS1413+135, open stars from B3 1504+377 at z=0.673 (Wiklind & Combes 1996b) and the filled triangle PKS1830-210 at 
z=0.886 (Wiklind & Combes 1996a). The dashed line is a one-to-one correspondance between the column densities and is not 
a fit to the data. 
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Fig. 11. Column density of HCN plotted versus column densitites of HCO + , HNC and CS. Designations are as in Fig. 10. 
The full-drawn line in N(HCN) vs. N(HNC) is a linear fit to the data, showing how the HCN / HNC ratio increases with 
decreasing HCN column density (see text for details). 



5.1.1. The extinction towards PKS1413+135 

The minimum column density of molecular hydrogen to- 
wards PKS1413+135, inferred directly from the CO(0->1) 
absorption, is 4 x 10 20 cm~ 2 (Sect. 4.1). Together with the 
estimated HI column density of 1 x 10 21 cm~ 2 (where we 
have assumed a covering factor of unity and a spin tem- 
perature of 100 K), the total column of hydrogen, N(H) = 



2N(H 2 )+N(HI), is - 2x 10 21 cm" 2 . This is more than 10 
times lower than the column inferred from the deficiency 
of low energy X-ray photons which indicates A v > 30 mag 
(Stocke et al. 1992). Is this discrepancy caused by an im- 
precise X-ray measurement, by the particular distribution 
of the obscuring gas relative to the radio core, or by an un- 
derestimate of the molecular and/or atomic column den- 
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sity? McHardy et al. (1994) analyzed ROSAT observations 
of PKS1413+135 in the 0.4-2.4keV range. Although they 
did not obtain enough data for a complete spectrum, the 
very low flux of low-energy X-ray photons measured by 
ROSAT is consistent with the deficiency observed by the 
Einstein satellite (Stocke et al. 1992). It is thus not likely 
that the X-ray results are erroneous. Also, the absence of 
any lines from both the narrow- and broad emission line 
regions, even when viewed in the near-IR (McHardy et al. 
1994; Stocke et al. 1992; Perlman et al. 1996), indicates 
that the extinction is very high. An obscuring dusty torus 
around the AGN can cause a large extinction. In CygA, 
for instance, the optical extinction, originating in a torus 
near the center, is estimated to be ^170 mag, or N(H) 
= 3.75 x 10 23 citT 2 (Ueno et al. 1994). Although 21cm HI 
absorption is detected (Conway & Blanco 1995), no molec- 
ular absorption is seen despite this high column density 
(Drinkwater et al. 1995). The non-detection of molecu- 
lar absorption could be due to either a lack of molecular 
gas as such, or that the gas is dense enough to render the 
excitation temperature very high, depopulating the lower 
rotational levels of existing molecules, or that the molec- 
ular rotational levels are radiatively coupled with the am- 
bient radiation field (e.g. Maloney et al. 1994). The pres- 
ence of a dense torus in PKS1413+135, similar to CygA, 
is, however, not likely in view of the absence of a near- 
infrared excess emission, caused by dust grains heated to 
high temperatures. This implies that the column density 
of the molecular gas in PKS1413+135, which should be 
situated at a relatively large distance from the AGN in 
order to explain the lack of near-infrared emission, and 
to be compatible with the very narrow line-width, is sev- 
erly underestimated, strongly suggesting the presence of a 
dense clumpy medium as emphasized in previous sections. 

5.1.2. Time variability 

Changes in the absorption lines towards background con- 
tinuum sources can be caused either by motions of small 
scale structure in the absorbing gas across the angular ex- 
tent of the radio core, or by spatial changes in the radio 
core itself. Such variations can give important informa- 
tion about both the size of the radio core and of the scale 
spectrum of the molecular interstellar medium. 

The radio source PKS1413+135 consists of a compact 
radio core and several emission components on a parsec- 
scale, while no radio emission is seen at scales of kilopar- 
secs (Perlman et al. 1994, 1996). The parsec-scale compo- 
nents all have a steep-spectrum and do not contribute to 
the flux at millimeter wavelengths. The core has a spectral 
index a — +1.7 between A18cm and A3cm (Perlman et al. 
1996). The core remains unresolved at an angular reso- 
lution of 2.3 mas, corresponding to ^7pc. It is, however, 
likely to be significantly smaller. Perlman et al. (1996) de- 
rive a lower limit to the size of the core from the variabil- 
ity time-scale of ~10/iarcsec, corresponding to ~0.03pc. 



The galaxy associated with PKS1413+135 is seen edge-on 
(McHardy et al. 1994) and we can therefore expect that 
the transverse velocity is equal to the rotational velocity. 
The latter is unknown but likely to be ~250kms _1 . This 
corresponds to a transverse shift of 50 AUyr -1 . The time 
scale for a significant change of the gas properties due to 
a shift of the obscuring molecular gas is therefore at least 
100 years, possibly much longer. The only possibility to 
explain much shorter variations is the existence of high 
velocity 25 000 kms -1 ) shocks propagating outwards 
from the radio core. 

In Table 4 we also give the continuum flux relative to 
that of May 1994 (which was the highest during the extent 
of these observations). The deepest absorption (largest 
ratio of the first and second component) occurs when 
the continuum flux is at maximum, suggesting that the 
changes are due to structural changes in the background 
radio source rather than transverse motion of small scale 
structure in the molecular gas. Knowing the extent of the 
changes of the radio core during an outburst would give us 
a limit to the size of the molecular gas structures. Future 
mm-VLBI observations may provide such data. Also, the 
change in component 1 is likely to be recurrent with the 
next outburst, if this takes place in the same region of the 
radio core. 

5.2. Molecular cloud properties 

The inferences to be made from the analysis of mul- 
ticomponent molecular gas along the line of sight to 
PKS1413+135 (Sect. 4.3 & 4.4) is that it is not possible to 
derive a correct column density without knowledge about 
the structure of the molecular ISM. However, if we can 
assume homogeneity for the physical and chemical condi- 
tions over the extent of the background source and along 
the line of sigth, we can quite confidently derive abun- 
dance ratiosLi. Furthermore, assuming similarities between 
the structure of the molecular ISM in PKS1413+135 and 
the Milky Way, we can compare our column densities with 
those derived for our Galaxy through similar observing 
methods (i.e. absorption line measurements). These as- 
sumptions are justified by the small amount of molecular 
gas that actually contributes to the observed absorption. 
The angular extent of the background continuum source is 
of the order 10-30 /larcsec, as implied by time variability 
(Perlman et al. 1996) and by inferences from mm-wave 
VLBI of a similar radio source (Lerner et al. 1993). The 
amount of molecular hydrogen probed by the line of sight 
to PKS1413+135 is thus only (0.5 - 2) x 10" 2 M Q , where 
we have used N(H 2 ) = 4 x 10 20 cm" 2 (cf. Sect. 4.1). The 
size-linewidth relation found for molecular clouds in the 
Milky Way (Larson 1981; Solomon et al. 1987) implies 
thay the molecular gas is extended ^lpc along the line 
of sight. Assuming that most of the mass is contained in 

3 Abundance is used in synonym with column density while 
discussing column density ratios. 
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Fig. 12. The HCN/HNC ratio versus the column density of HCO+, HCN and HNC. The trend of decreasing HCN/HNC 
abundance ratio with increasing column density of HCN indicated in Fig. 11 is clearly visible. Designations are as in Fig. 10. 



small clumps with n(-ff 2 ) — 10 4 cm~ 3 (Sect. 4.4), the to- 
tal H 2 mass is /y(0.4 — 4)M Q , where /y is the volume 
filling factor of the dense clumps, jy is likely to be <<1. 
In either case we are dealing with a very small portion 
of the molecular ISM in the galaxy at z=0.247. In Ta- 
ble 5 we present abundance ratios for PKS1413+135. The 
column densities are taken from Table 2. The upper and 
lower limits represents la limits. Below we discuss in some 
more detail various results that can be derived from these 
abundance ratios. 

5.2.1. Molecular oxygen O2 

One of the more interesting questions in interstellar chem- 
istry is the abundance of 2 , which is supposed to bind 
most of the free oxygen atoms and to be one of the most 
important coolants (e.g. Goldsmith & Langer 1978). This 
molecule has yet to be observed in the ISM. The redshift 
of PKS1413+135 shifts the main N(J)= 1 -> li line of 
O2, which is usually not observable from the ground, to an 
easily accesible window (95.25 GHz) However, the rather 
low average opacity means that our limit to 02/CO< 
0.3 ± 0.09 is 20 times higher than the upper limit derived 
towards B0218+357 at z=0.68 (Combes & Wiklind 1995). 
Nevertheless, this limit is of interest since it concerns the 
ground state and thus is sensitive to excitationally very 
cold gas. The transitions observed in B0218+357 were the 
N(J)= 1 2 -> 3 2 and the N(J)= U 3 2 lines. The li 
level can not be directly populated from the ground state 
through collisions, but has to be populated through col- 
lisional excitation from N(J)= lo to 1 2 and 3 2 (Bergman 
1995). The latter energy level corresponds to a tempera- 
ture ~23K. Radiative transitions from the ground state 
to the N(J)= li level is allowed and probed by the obser- 



vations presented here. This transition has now also been 
observed in B0218+357, giving an improved upper limit 
to the 2 /CO ratio for this system (Combes & Wiklind 
1997). 
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5.2.2. Isotopic ratios 

The isotopic abundace ratios for CO and HCO + presented 
in Table 5 suggests relatively low abundances of the 13 C 
isotopic variants. The large formal error for 12 CO/ 13 CO 
comes from the assumed range of excitation temperatures 
(cf. Table 2). Since the abundance ratio is largely inde- 
pendent of T x , the true uncertainty of this ratio is much 
smaller. In the ISM of the Milky Way the 12 C/ 13 C isotope 
ratio varies from —20 in the Galactic center region to ~70 
in the local ISM (Wilson & Matteucci 1992). The correla- 
tion between the 12 C/ 13 C ratio and that of the substituted 
isotopic molecules depends on chemical fractionation (en- 
hancing the isotopic variants) and selective photoinization 
(decreasing the isotopic species due to less self-shielding) . 
In Galactic molecular clouds as well as in nearby galax- 
ies, the final result is an increase in the 12 CO/ 13 CO ratio 
(and likewise for HCO + ) compared to the actual 12 C/ 13 C 
ratio. The high lower limits found in PKS1413+135 there- 
fore suggests a low 13 C abundance relative to 12 C. 

Whereas 12 C is produced in both low and high mass 
stars, as a primary product of hydrostatic burning, 13 C 
is only produced through incomplete proton burning in 
the red giant stage of low and intermediate mass stars 
(e.g. Wilson & Matteucci and references therein). A low 
13 C abundance may therefore indicate a young chem- 
istry, where the low and intermediate mass stars have not 
yet had time to reach the red giant stage. 14 N is pro- 
duced and ejected into the ISM in much the same way 
as 13 C. The high lower limit HCO+/N 2 H+ > 41 found 
in PKS1413+135 is therefore in aggrement with the inter- 
pretation that the low and intermediate mass stars have 
not yet reached their red giant stage. In PKS1830-211 at 
z=0.89, we found a HCO+/N 2 H+ ratio of -1.4 (Wiklind 
& Combes 1996a). In this molecular absorption line sys- 
tem we also detect several 13 C isotopic species. Although 
it is not entirely clear whether this is due to an extremely 
high opacity, PKS1830-211 seems to be more chemically 
evolved than PKS1413+135, despite the much higher red- 
shift. 

5.2.3. Dark cloud chemistry 

The equilibrium abundances in dark clouds (i.e. molecular 
gas where the main ionization source is cosmic rays) can 
be divided into a High and Low Ionization Phase (HIP and 
LIP) (Flower et al. 1994; Le Bourlot et al. 1995; Gerin et 
al. 1997). The LIP is characterized by a chemistry driven 
by ion-molecule reactions where proton transfer involving 
Hjj~ plays a major role. The C/CO ratio is low and most of 
the gas phase carbon is locked up in CO. This phase has a 
high abundance of 2 and molecular ions such as HCO + 
and N 2 H+. Other carbon bearing molecules, such as CN, 
have low abundances. The HIP has a chemistry driven by 
charge transfer reactions involving H + and is characterized 
by a large C/CO ratio. This leads to low abundances of 



O2 and molecular ions. Abundances of HCN and HNC 
are only marginally affected by the ionization state of the 
gas. The abundance ratios presented in Table 5 mainly 
uses HCN as a 'reference'. 

Are our abundance ratios for PKS1413+135 consistent 
with either of these two phases for dark cloud chemistry? 
Comparing our abundance ratios with the model results 
of Le Bourlot et al. (1995) we find that our O2/CO ra- 
tio docs not constrain the ionization state of the gas, al- 
though our limit is of the same order as the expected 
O2/CO ratio for the low ionization phase. Our limit to 
the O2/HCN ratio is only compatible with a high ion- 
ization phase, while our abundance ratios for HCN/CO, 
HCN/HNC and the lower limit to HCN/CS only are com- 
patible with a low ionization phase. Hence, the molecular 
absorption line data for PKS1413+135 does not present 
a clear-cut case for either of the two ionization phases. 
Also, the HCO+/CO ratio is too large, by more than an 
order of magnitude, to be compatible with either the high 
or low ionization phase. However, as discussed by Hoger- 
heijde et al. (1995), turbulence might influence the chem- 
istry, producing enhanced amounts of HCO + . If this is 
the case for the high abundance of HCO + this would indi- 
cate a diffuse gas, since the proposed formation pathway 
is quenched at n(i?2) ^ 10 3 cm~ 3 . Dark cloud chemistry 
shows bi-stability and can switch from one phase to the 
other on a short time scale (Le Bourlot et al. 1992, 1995), 
making it possible to have a mixture of ionization phases 
over small spatial scales. Improved limits on the O2 and 
CS abundances may enable us to put limits to the amount 
of gas along the line of sight to PKS1413+135 that can 
exist in HIP and LIP states. 

5.2.4. Comparison of column densitites 

In Fig 10 we plot the column density of HCO + versus the 
column densities of HCN, HNC and CS. The plot contains 
data from Galactic absorption measurements in low den- 
sity gas (Lucas & Liszt 1993, 1994, 1996) and in higher 
density gas seen towards SgrB2 (Greaves & Nyman 1996), 
as well as our recent results from absorption line measure- 
ments towards Cen A (Wiklind & Combes 1997). In addi- 
tion, we include results from two additional high-z molec- 
ular absorption line systems 1504+377 at z=0.67 (Wiklind 
& Combes 1996b) and PKS1830-211 at z=0.89 (Wiklind 
& Combes 1996a). The dotted line is a one-to-one corre- 
spondance between the column densities and not a fit to 
the data. In Fig. 1 1 we plot the column densities of HNC 
and CS versus HCN in the same way as in Fig. 10. Despite 
the presence of a considerable scatter, the most striking 
impression from the figures is the remarkable correlation, 
over more than 3 orders of magnitude, in column density. 
This suggests that the chemistry is similar for these simple 
molecules in environments characterized by very different 
column densities. The high-z systems do not show any pe- 
culiarites compared to the local values, suggesting that the 
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molecular ISM offers similar conditions for star formation 
at earlier epochs as it does in the present one. 

The only trend in abundance ratios is found between 
HCN and HNC (Fig. 11), where the abundance of HNC 
seems to decreases relative to HCN with decreasing HCN 
abundance. The solid line is a least-square fit to the HNC 
vs. HCN abundances: 

N(HNC) = 2.3 x 1(T 3 [N(HCN)] 1A6 . (10) 

In Fig. 12 we plot the HCN/HNC abundance ratio as a 
function of the HCO + , HCN and HNC column densities. 
The HCN/HNC ratio decreases with increasing column 
density in all three cases. Isomers like HCN and HNC 
have similar chemistry only differing in a few key reactions 
which determine the abundance ratio. Through chemical 
modelling of the HCN and HNC abundances Schilke et 
al. (1992) found that the HCN/HNC ratio increases with 
increasing kinetic temperature (cf. Irvine et al. 1987) as 
well as with increasing molecular hydrogen density, nn 2 . 
The reason for this is primarely the key reactions 

HNC + O — > NH +CO 
HNC + C — > HCN + H 

which are effective at relatively high temperatures and 
preferentially destroy HNC relative to HCN. A large col- 
umn density does not necessarily imply a large volume 
density of H2 and, since the abundance of HCN is less 
dependent on Tk for low H 2 abundances (Schilke et al. 
1992), the correlations between the HCN/HNC and col- 
umn densitites seen in Fig. 12 is likely to be a temperature 
effect, affecting primarily the HNC abundance. The low- 
est column densities being associated with the kinetically 
warmest gas. However, the dependence of the HCN/HNC 
abundance ratio on both temperature and density makes 
it less valuable as a temperature probe then previously 
thought. 

The high-z systems all tend to have HCN/HNC ratios 
~2, among the lowest of the line-of-sights presented in 
Figs. 10-12. It is unclear whether this is a selection bias 
towards both kinetically and excitationally cold gas, or 
an effect of low metallicity, making the HNC depletion 
reactions less effective. 

5.3. Invariance of physical constants 

A potential diagnostic application of molecular absorp- 
tion lines at high redshift is to check the invariance of 
fundamental physical parameters, (cf. Thompson 1975). 
The energy difference between two adjacent rotational 
levels in a molecule consisting of two atoms is propor- 
tional to (^r 2 ) -1 , where r is the bond length and 1S 
the reduced mass of the atoms. A vibrational transition 
of the same molecule has, to a first order approximation, 
a ^i 1 / 2 dependence on the reduced mass. Comparison of 



Table 6. Redshifts and variations of molecular mass 
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a) Az/(1 + z) = (zhi - z mo i)/(l + z HI ) 

HI references: PKS1413+135 (Carilli et al. 1992), B0218+357 (Carilli 
et al. 1993), B3 1504+377 (Carilli et al. 1996). 

Molecular line references: PKS1413+135 (this paper), B0218+357 
(Combes & Wiklind 1995), B3 1507+377 (Wiklind & Combes 1996b), 
PKS1830-211 (Wiklind & Combes 1996a,1997). 



a vibrational-rotational spectrum from an earlier epoch 
with a present one can therefore be used to set limits on 
the changes of the proton and neutron masses. This has 
been done for molecular hydrogen seen in absorption in 
a damped Lyman-a system at z=2.81 (Foltz et al. 1988, 
Varshalovich & Potekhin 1996), giving a limit of 2 x 10 -4 
to Am p /m p . For pure rotational transitions one can com- 
pare the line frequencies (or similarly, the measured red- 
shifts) with the frequencies of electronic transitions. This 
gives a measure of the invariance of the electron to proton 
mass. It is also possible to compare the 21 cm HI line with 
a molecular rotational line. The formers energy difference 
to first order depending on m 2 , again giving a measure 
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Fig. 13. A comparison of the CO(0— »1) and 21 cm HI absorp- 
tion towards PKS1413+135. In both cases the optical depth is 
plotted. Notice that Fig. 3 in Wiklind & Combes 1994 suffers 
from an erroneous velocity scale, causing the CO and HI lines 
to be separated by ~llkms _1 . 



of the invariance of the m e /m p ratio (see also Tubbs & 
Wolfe 1980 for a comparison of the 21 cm line with optical 
resonance transitions) . 

In an Appendix we express the energy levels of molec- 
ular rotational lines in fundamental physical parameters. 
It is shown that for a differential measurement, i.e. mea- 
suring the redshifts of two different molecular species, 
Az/(l + z) is directly proportional to A^i/fi. Strictly 
speaking, this is only valid when the molecule is approxi- 
mated as a rigid rotator, but this is a good approximation 
for low lying J states. In Table 6 we list the observed red- 
shifts of the molecular absorption lines in PKS1413+135, 
together with some lines from B 0218+357, B3 1504+377 
and PKS1830-211. We also list the redshifts derived from 
21cm HI absorption when available. The redshifts have 
been obtained using Gaussian fits. We can compare both 
the redshifts derived for different molecules and the red- 
shifts of the HI and molecular lines. If changes in the elec- 
tron mass are negliblc, the Az/(l + z) measures A/i//i, 
to first order, in both cases. The values given in Table 6 
are the difference between the 21 cm HI and molecular 
redshifts. In Figure 13 we show the 21cm HI absorption 
(Carilli ct al. 1992) together with our CO(0->1) absorp- 
tion towards PKS1413+135. It is quite clear that the HI 
spectra limit the accuracy with which Az can be mea- 
sured. Except for B3 1504+377, all the values in Table 6 
are compatible with Az = 0. For the 21cm HI absorp- 
tion obtained towards B3 1504+377 (Carilli et al. 1997), 
the discrepancy between HI and molecular velocities is 
of the order of SOkms" 1 , i.e. Az « 10~ 4 , an order of 
magnitude larger than the HI accuracy. The HI result has 
been confirmed to an accuracy better than lkms -1 with 



the Westerbork Synthesis Radio Telescope (Carilli et al 
1997, preprint). Hence we reach here the basic limitation 
of the method, consisting of observing different atoms or 
molecules. It cannot be expected that all lines peak at 
the same velocity, since their local association is not per- 
fect, and their excitation might be different. Supporting 
this, note that the molecular component in absorption in 
B3 1504+377 at z = 0.67150 has no HI counterpart. 

To obtain more statistics, atomic and molecular ab- 
sorptions in the Milky Way must be considered. Liszt & 
Lucas (1997) compared the absorption profiles of several 
different molecular lines with that of 21cm HI and found 
that the velocity of the peak opacity often differs between 
the molecular and atomic lines. This shows that the ab- 
sorptions do not arise in the same gas and any comparison 
is limited by the possible velocity difference of the two ISM 
components where the absorptions originates. 

It is also possible to use different transitions from the 
same molecule to estimate A/if/i; While the rotational 
constant is proportional to yU _1 , the centrifugal distortion 
coefficient is proportional to ji~ 2 (Towncs & Schawlow 
1975). In the Appendix we show that we actually have 
a dependence on the invariance of the electron mass and 
the fine structure constant in this case. Comparing the 
J=1^2 and J=2^3 lines of HCO+ in PKS1413+135 we 
get at the la level: Az/(l + z) = (4.0 + 3.4) x 10~ 7 . Simi- 
larly for HCN: Az/{l + z) = (5.6 + 8.4) x 10" 7 . Comparing 
CO(0->1) and HCO+(1^2) we get: (-1.1 ± 0.3) x lO" 6 . 
This last result is marginally different from zero at a 3a 
level. However, even when comparing molecular transi- 
tions one has to face the possibility that the lines do not 
come from the same location in phase space. HCO + , for 
instance, has been shown to exist in regions where the CO 
molecular is no longer self-shielded (Lucas & Liszt 1996), 
and its absorption profile can therefore be weighted with 
velocity components not seen in the CO line. A further 
complication is that HCO + is an ion and susceptible to 
magnetic fields present in the ISM. This can give it a drift 
velocity different from neutral molecular species. It is pos- 
sible to have absorption lines dominated by different ve- 
locity components even when considering transitions from 
the same molecular species. This is due to the nonlinear 
dependence of the opacity to the excitation temperature 
(Sect. 4.2). Different transitions can therefore be weigthed 
differently by gas along the line of sight. 

To summarize, measurements of rcdshift differences 
between molecular absorption lines and 21 cm HI are likely 
to be dominated by velocity differences between the two 
gas components along the line of sight. This is a severe 
concern when comparing different molecular species and 
even different transitions of the same molecule. Taken at 
face value, the HI accuracy limits the nucleon mass varia- 
tion to typically Afi/fi <J 10 -5 at a 3a level. In one source 
(B3 1504+377) a discrepancy as large as 10 -4 has been 
observed. 
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6. Summary 

We present a multi-line study of an absorbing molecu- 
lar cloud in the edge-on galaxy hosting the radio-source 
PKS1413+135, at z=0.24671. We have shown that the to- 
tal column densities obtained are lower limits, since the 
absorption essentially traces the diffuse extended molec- 
ular gas, and is only weakly weighted towards the dense 
phase which provides most of the mass. The size of the ab- 
sorbing structure along the line of sight must be of parsec 
scale, given its narrow total width. This situation is favor- 
able to detect time variations in the absorption profile, if 
there exists knots in the radio source moving close to the 
velocity of light. Variations on time-scales of a month then 
correspond to structures of > 10 3 AU. We report possible 
time variation, affecting the relative ratio of two absorb- 
ing components, and therefore the shape of the profile. 
This is to a large extent independent of the continuum 
level and calibration uncertainties. Such time-variations 
are not surprising, in view of the small scale structure 
that has been detected in the ISM of the Milky Way. 

Because of the high spectral precision of the millime- 
ter techniques, it has been proposed that observation of 
molecules at high redshift could constrain variations of 
fundamental physical constants, such as the mass of the 
molecules. However, different molecules or different kind 
of lines (rotational, vibrational, etc..) should then be com- 
pared. Concerning the comparison between HI and molec- 
ular lines, it was shown that the HI accuracy limits the 
nucleon mass variation to A/x/^u ^ 10~ 5 at 3a in many 
cases, but a discrepancy as large as 10 -4 has been ob- 
served in one of the sources. Excitation problems in the 
not well-known multiphase structure of the medium pre- 
vents reaching conclusive results. 
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7. Appendix 

Here we present in more detail the use of molecular rotational 
lines as a probe of the invariance of physical constants over 
large spatial and temporal scales. In particular, we show that 
there is a limit to the accuracy of molecular rotational lines, 
coming from our incomplete solution of the Schrodinger wave- 
function for a complex system such as a molecule. 

7.1. A diatomic rigid rotator 

As a first approximation we will look at a diatomic molecule 
without considering the centrifugal distortion and the zero- 
point vibrational energy. 

The energy of a rotational level J is 



Erot = -£yJ(J+l) 



(11) 



where J is the angular momentum quantum number and I 
is the molecular moment of inertia about axes perpendicular 
to the internuclear axis (cf. Townes & Schawlow 1975). For a 
diatomic molecule the moment of inertia can be expressed as 
pjr\ , where n is the reduced mass of the two atoms and r c is the 
equilibrium intranuclear distance of the two atoms. To within 
a factor of 2-3, we can set ft~ao, the Bohr radius. Expressed 
in fundamental parameters 4 !, the energy of a rotational level J 
is 



Ei, 



/' 



-J(J+ 1) 



(12) 



where m c is the electron mass and a the fine structure constant. 
The frequency of a rotational transition following the selection 
rule AJ = ±1 is 



2^ / r 1 \ I c 



(13) 



This gives the correct frequency of a transition to within a fac- 
tor of a few. The discrepancy arises from the approximation 
of the effective internuclear separation as the Bohr radius, but 
does not affect the analysis of the fundamental physical pa- 
rameters. A shift in the line frequency due to changes in the 



4 These are usually called 'constants', but since we are testing 
their invariance, we will refer to them as 'parameters'. 



reduced mass fi, the electron mass m e and the fine structure 
constant a is therefore 



2 

+ 4 f 



(14) 



If the redshift of an object causing the observed transition is 
known apriori, we can directly relate the shift in frequency or, 
similarly, in velocity Av = cAvjv to combined changes in the 
three parameters (i, m c and a. In order to measure a shift of 
1 kms - 1 , we must have Afx/fj, ~3x 10~ 6 (if Am c = Aa = 0). 
However, any determination of a redshift involves a measure- 
ment of line radiation, with the possibility of being affected by 
changing physical parameters (not necessarily the same as for 
molecular rotational transitions). Hence, we must make a 'dif- 
ferential' measurement in order to test their invariance. In the 
case of a rigid diatomic molecule, this is done by comparing two 
different species. What we measure is two ratios of frequencies; 
the 'laboratory' (called lab) and the observed frequency (called 
obs): 



flab 
fobs 



1 + 2. (15) 
For two different molecules a and b, we can now write: 



"lab 






"obs 


b 


Az 


"lab 




1 + Za 


"obs 


a 





(16) 



Since we are comparing two molecular species at two different 
epochs, the electron mass and the fine structure constant will 
drop out and we have 



Az 

1 +Za 



ARy 



A-R u = I — 



Ru 



(17) 



(18) 



(19) 



In Eqs. 20-22 we use as a subscript for mass ratios obtained in 
the laboratory (z=0) and z for mass ratios at a distant epoch. 
Comparison of molecular rotational transitions are thus sensi- 
tive to the strong force only, affecting the proton and neutron 
masses. For practical purposes, the largest effect is obtained 
when J? M0 is small, i.e. >> Mb- 

7.2. Non-rigid diatomic molecules 

The next step in the complexity of a diatomic molecule is 
the centifugal distortion and the zero-point vibrational energy. 
As the rotational quantum number J increases the molecule 
rotates faster and the internuclear separation increases. This 
leads to a larger moment of inertia and a decrease in the energy 
compared with the rigid rotator. The zero-point vibrational 
energy introduces a periodic oscillation in the internuclear sep- 
aration, affecting the moment of inertia and the energy levels. 

A separate treatment of the atoms and the electrons (the 
Born-Oppenheimer approximation) is necessary in order to de- 
scribe the total potential in such a way that the Schrodinger 
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wave function can be solved. Ultimately, as we will discuss 
in the next subsection, the Born-Oppenheimer approximation 
puts a limit to the accuracy with which molecular rotational 
transitions can be used to test the invariance of physical pa- 
rameters. 

Dunham (1932, Phys Rev 41, 721) gave the vibrational 
and rotational energy levels for any potential which can be 
expanded in powers of (r — r e ) (see also Townes & Schawlow 
1975). However, we will use the simpler Morse potential in 
order to emphasize the dependence on fundamental physical 
constants. The Morse potential is written 



U(r) = D(l - exp(-a(r - r c )) 2 



(20) 



where D is the dissociation energy of the molecule, r c is the 
equilibrium intranuclear separation and a is a constant. The 
frequency of a rotational transition (v = 0) becomes 



2B e (J + 1) - 4D e (J + 1) - a e (J + 1) 



n 



47Yl e 

8ttVb| 
a 2 D 

3c 3 



2h 



(21) 

(22) 

(23) 

(24) 
(25) 



If we for simplicity disregard the zero-point vibrational energy, 
we can write the line frequency as 



7V 2 C 6 \ 

a 2 Dh :i J 
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(26) 



The first term is identical to that obtained for the rigid rotator 
and the second term is the effect of centrifugal stretching (low- 
ering the frequencies). Due to this term, the dependence on m e 
and a at z — and at z does not cancel when comparing the 
redshifts for two molecules a and b, as was the case for the rigid 
rotator. Hence, Az/(1 + z) depends on all three parameters; 
fj,, m c and a. For low J-values, however, the dependence on 
the electron mass and the fine structure constant is small, and 
the rigid rotator represents a good first order approximation. 
The nonlinear dependence on J in the centrifugal distortion 
means that it is possible to compare the redshifts derived from 
different rotational transitions of the same molecule. 



There are two main uncertainties which affect a test of 
the invariance of fundamental physical parameters. The first 
is anharmonicity of the potential. This is particularly true for 
the Morse potential used here. The situation is improved some- 
what by using the method of Dunham (cf. Townes & Schawlow 
1975), but the method involves parameters which can be cal- 
culated from spectroscopically observable quantities and they 
ultimately put a limit to how accurate the ratio of two re- 
duced masses can be calculated to ~ 10~ 7 (Townes & Schawlow 
1975). The second uncertainty involves the Born-Oppenheimer 
approximation, where the electronic and atomic contributions 
to the total potential are treated separately. Although electrons 
can be considered to be more or less spherically distributed 
around their respective nuclei, and thus not contributing to 
the moment of inertia (the so called 'slip effect'), electrons in 
the valence shell will give a contribution to the moment of 
inertia of ~ n c m c r 2 , where n c is the number of electrons in 
the valence shell, m c is the electron mass and r is an average 
radius from the nucleus with which the valence shell is associ- 
ated. This contribution to the moment of inertia is very small 
indeed, but when deriving mass ratios it limits the accuracy 
to ~ 10~ 6 (smaller for heavy molecules and larger for lighter 
ones) . An effect which is closely related to the valence electrons 
is so called L-uncoupling. Rotation tends to excite the valence 
electrons from their ground 1 E state with zero angular momen- 
tum, to excited 1 TL states, with unit angular momentum. This 
in itself affects the rotational energy, but a larger effect comes 
from the fact that the n state also produces a large magnetic 
field at the positions of the nuclei. This induces magnetic hy- 
perfine interaction, again changing the rotational energy levels. 
This L-uncoupling can be measured in the laboratory and is 
found to introduce uncertainties of ~ 10~ 6 in estimated mass 
ratios (Townes & Schawlow 1975). 

Hence, the uncertainites associated with estimates of mass 
ratios of molecules is of the order 10~ 6 . This value is of the same 
order as the A/i/n necessary to produce a shift of ~lkms _1 
relative to the 'true' redshift of a molecular rotational line. 
Smaller shifts than this will therefore not be useful in testing 
the invariance of physical parameters. 



7.3. Limitations to the accuracy 

As discussed in Sect. 5.4, the accuracy of using comparative 
measures of redshifts derived from either different emission 
mechanisms (optical, 21 cm hyperfine line, molecular rotational 
lines, etc) is limited by velocity gradients in the emitting 
medium and the fact that different gas masses contribute dif- 
ferently to the various emission/absorption processes. This is 
a grave concern also when using molecular rotational transi- 
tions, even when the same molecular species is used, but with 
different J levels. However, even if this was not the case, the 
accuracy is limited by our inexact quantum mechanical de- 
scription of a molecule. 
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